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To probe naked spacetime singularities with waves rather than with particles we study 
the well-posedness of initial value problems for test scalar fields with finite energy so that 
the natural function space of initial data is the Sobolev space. In the case of static and 
conformally static spacetimes we examine the essential self-adjointness of the time trans- 
lation operator in the wave equation defined in the Hilbert space. For some spacetimes 
the classical singularity becomes regular if probed with waves while stronger classical sin- 
gularities remain singular. If the spacetime is regular when probed with waves we may 
say that the spacetime is "globally hyperbolic." 
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1 Introduction 



In general relativity a singular spacetime is defined by the geodesic incompleteness Q. 
However, sometimes such a definition gives a very weak singularity which seems almost 
harmless from physical point of view. For example, a spacetime from which a single point 
is taken out is a singular spacetime because there is a geodesic curve which terminates 
at the point outside of the spacetime with a finite affine time. A stronger 'physical 
singularity' appears for example at the center of a spherically symmetric black hole, where 
the curvature scalar diverges and therefore the resultant infinite tidal force will tear off 
any physical object. The classification of singularities is yet under way but has not been 
completed 0. 

The standard definition of a spacetime singularity is physically based on a probe with 
classical point particles. In this paper we shall discuss a wave probe of timelike singulari- 
ties which was initiated by Wald pi and later developed by Horowitz and Marolf M . The 
idea of the probe with waves rather than with classical particles is motivated by quantum 
field theory because everything should be described by quantum fields. The wave may 
propagate through the would-be singularity with a definite and unique way. For example, 
in the case of hydrogen atom the wave function is finite at the origin, which is a classical 
singularity. It is known that if the space is geodesically complete the Laplacian operator 
has a self-adjoint extension and the extension is unique so that the wave propagation is 
well defined. Converse is not always true. If the geodesic completeness is replaced by the 
well-posedness of initial value problems for test fields the concept of the global hyperbol- 
icity and therefore the cosmic censorship f] should be drastically changed as Clarke [pfl has 
advocated. 

We shall be concerned with a wave propagation dictated by the Klein-Gordon equation 
in a curved spacetime with timelike singularities. Only for illustration in the introduction 
we use the simplest case; the Klein-Gordon equation in (1 + l)-dimensional Minkowski 
spacetime, (— df + d 2 )f = defined in a suitable region of the spacetime. (For a general 
case, see the following sections.) 

For the initial value problem we introduce the following norm on a function space on 
each t = const, hypersurface: 

imh= (jJ dx \f\ 2 + lJ dx 

where q 2 is a positive constant. We call the function space W = {/|||/||<oo} the Sobolev 

3 Here we refer to the "physical formulation" of the strong cosmic censorship rather than the "precise 
formulation" in Wald's book If our wave approach can be extended to the initial value problem of 
the Einstein equations the notion of the cosmic censorship will substantially change. 



df 
dx 
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space or H 1 . The Sobolev norm has been used in the standard formulation of well-posed 
initial value problems in a general globally hyperbolic spacetime 0. We note that in 
general the well-posedness of an initial value problem requires continuous dependence 
of solutions on initial data, besides the existence and the uniqueness of solutions ||. 
However, the main issue we will address in this paper is to see the uniqueness of solutions 
of a wave equation in a non-globally hyperbolic spacetime, so hereafter we say that the 
initial value problem is well-posed when the wave propagation is uniquely determined in 
the whole spacetime. [] 

It is known that the norm is bounded above by the field energy so that the finiteness 
of the energy implies the finiteness of the norm 0. Since we cannot afford to prepare an 
infinite energy field configuration as initial data, the function space is naturally limited 
by the condition: [] 

ll/li<oo. (2) 
The corresponding natural inner product in our Hilbert space is defined as 

(f, 9) :=£j^ rg + y^£ d £ (3) 

so that Il/H 2 = (/,/). 

We shall confine ourselves mainly to the case of timelike singularities in static or 
conformally static spacetimes so that the wave equation becomes of the form, 

d 2 A = -A(f>, (4) 

where A is an operator which contains spatial coordinates and spatial derivatives only. In 
this case the well-posedness of the initial value problem is translated into the essential self- 
adjointness of the operator A because of the spectral theorem ||. Namely, we prepare 
a smooth and nice initial data at some spatial hypersurface by choosing the Sobolev 
space H 1 as the Hilbert space. The wave will propagate and eventually hit the timelike 
singularity and will be scattered off in someway. The point of the essential self-adjointness 
is that any unwanted singular modes which are not contained in the domain of the initial 
data will not appear after scattering so that the initial value problem is well-posed with 



4 To show the existence of solutions and to establish an appropriate continuous relation between initial 
data and solutions, Sobolev norms containing higher order derivatives are chosen to define a topology on 
the space of initial data. However, to prove the uniqueness of solutions of second order linear hyperbolic 
equations, it is sufficient to adopt H 1 as our Sobolev space which is larger than H m with norms 
containing m(> l)th order derivatives. Our results in this paper hold also for H m ^ >1 \ 

5 The difference between ours and that of reference is in the definition of the norm and therefore 
of the Hilbert space. In the case of quantum mechanics the natural Hilbert space is the linear function 
space with the square integrability, L 2 , because of the probabilistic interpretation of the wave function. 
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no arbitrariness in the choice of the boundary conditions and the prediction is unique. In 
such a case we say that the spacetime is "wave-regular." 

We would like to emphasize the relevance of the present work to quantum field theory 
in curved spacetime. It will be natural to expand a quantum field in terms of the nor- 
mal modes which belong to the Sobolev space rather than L 2 . We assign the coefficients 
of the mode expansion as annihilation and creation operators. The quantum states are 
constructed by applying the creation operators to the vacuum state which is defined by 
the condition that the vacuum is annihilated by all the annihilation operators. f\ This 
construction implies that if the initial value problem is well-posed the vacuum expec- 
tation value of the energy momentum tensor should be well-behaved near the would-be 
singularities so that the field energy is finite. 

The organization of the rest of the paper is as follows. In Sec. |2] we propose a natural 
choice of the function space in which the initial value problem is explored (Subsec. |2.1|) 
and we recapitulate the criterion of essential self-adjointness of operators in the Hilbert 
space (Subsec. |27|). In Sec. |]we demonstrate how we can probe singularities with waves 
in Minkowski spacetime with a single point removed and give an intuitive justification of 
the choice of the Sobolev space as the Hilbert space. Section f| supplies several examples of 
static spacetimes with timelike singularities. We explicitly show that many of the classical 
singularities become wave-regular, while a single example is wave-singular. Section |5| is 
the extensions of the discussion of the previous sections to scalar fields with general non- 
minimal coupling and to conformally static spacetimes. In Sec. [5] we discuss how to 
characterize wave-singular naked singularities in our approach and propose a notion of 
hair of naked singularities. Section [7| is devoted to summary and discussion. In Appendix 
some mathematical materials on the essentially self-adjointness are given for the sake of 
reader's convenience. 

2 The function space of initial data 

2.1 Finite energy field configuration 

We consider an (n + 2)-dimensional static spacetime of the metric form 



ds 2 = -V 2 dt 2 + h i jdx i dx i , 



(5) 



with a timelike Killing vector field £ M = (dt)^. 

We choose a function space on each t = const, hypersurface E as 




(6) 



6 The constructed quantum states belong to L 2 class in the Fock space. 
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with the Sobolev norm ||/|| being given by 

||/|| 2 :=^ JjEV- 1 /*/ + ^ JjZVhVDifDjf, (7) 

where q 2 is a positive constant and Di denotes the covariant derivative with respect to 
the induced metric hij on E. Here cffi = d n+1 x\/h is the natural volume element on S. 
The norm is bounded above by a positive constant times the energy integral E, 

\\f\\ 2 < const, x E, (8) 

where 

E:= J^dEn^TM (9) 
with being the unit normal to S. Here the energy momentum tensor is given by 

TM ■= ~ (V»f*VJ + VJ*VJ) - l -g, v (V7*VJ + m 2 r/) • (10) 

For = V~ 1 (dt) fl the energy E is expressed by 

E = \jj?. (V^dJ^J + m 2 Vrf) + l - jf dEVWDifDrf, (11) 

which motivated us to choose the norm given by Eq. (|^). The finiteness of the norm, 
ll/H < oo, is required because we can prepare only a finite energy configuration of the 
field. This leads us to the Sobolev space as the function space Ti on S. The energy 
E is conserved because the energy momentum tensor T^ v satisfies the conservation law: 
V^T 7 ^ = and £ M satisfies the Killing equation: C^ds 2 = 0. Then the inner product is 
naturally defined by 

(Lg):=^JjXV- 1 f*g + yjZVtiW i rD j g. (12) 

We will consider the massless case only because it is known that the initial value 
problem is well-posed for m ^ if it is for m = [|]. 

2.2 Uniqueness of the time translation operator 

Let us briefly recapitulate the mathematics on the essential self-adjointness of a linear 
operator A on the Hilbert space Ti. For precise definitions see Appendix, in which we 
collect relevant mathematical materials. 

7 Of course, the converse is not necessarily true. That is, ||/|| 2 < oo does not mean that the energy is 
finite in general. However, in our present analysis of the spacetime with a timelike (conformal) Killing 
vector, the Sobolev space implies the finiteness of the field energy. 
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The wave equation of a massless test scalar field, □</> = 0, reduces to 

d 2 (f) = -A*, (13) 

where A := —VD t VD i is a positive symmetric operator on 7i if the domain of A is 
suitably chosen, e.g., C^°(S), a set of smooth functions with compact support on E, so 
that it is dense in TC. In other words we see by a simple computation, 

(Af, g) = (f, Ag) + dS l {(A + q 2 )f Vd i9 - Vd t J*(A + q 2 )g} (14) 

so that A is symmetric if /, g £ C^°(S) and therefore the surface term above vanishes. In 
most cases this choice of the domain is not very restrictive. 

The domain of A can be further extended by relaxing the boundary condition so that 
the extended domain coincides with the domain of its adjoint operator. The extended 
operator in this manner is said to be self-adjoint and its eigenvalues are real and positive. 
Then, for each self-adjoint extension Ae, the time evolution of the field is uniquely given 
by i 

<p(t) = cos(A^ 2 t)(f)(0) + A~ 1/2 sin(^ /2 t)0(O), (15) 

with (f)(0), (f)(0) G T>(Ae) being any initial data. In this sense the self-adjoint extension 
Ae is a time translation operator. 

If there are many possibilities of the self-adjoint extensions, we have to choose one of 
them by imposing a particular boundary condition, which is normally imposed by some 
physical requirement. In the case of naked singularities we do not have any criterion to 
choose the boundary condition. Therefore, if the self-adjoint extension is unique, there 
remains no ambiguity in the choice of the boundary conditions. A symmetric operator A 
which has a unique self-adjoint extension is called essentially self-adjoint. 

The well-posedness of the initial value problem of Eq. flT3"D is now turned into the 
essential self-adjointness of the operator A, which can be tested by considering solutions 
of the equations 

A*ip = ±iip, (16) 
and showing that such solutions do not belong to our Hilbert space ||. 

3 Space with a single point removed 
3.1 Solution of the wave equation 

Let us consider a rather artificial model of a timelike singularity which can be fully 
analyzed. Namely, we consider a spacetime which is locally flat but with a single spatial 
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Figure 1: A conformal diagram of a spacetime with a timelike singularity at the center. 



point removed so that the spacetime has a timelike singularity as illustrated in Fig. [I] and 
the topology is (R 3 — {0}) x R. Our problem in this section is to see the well-posedness 
of the initial value problem of the Klein-Gordon equation, 

-<9 t 2 + A0 = O, (17) 

in this spacetime, which hopefully enhances our understanding of the wave probe for more 
general timelike singularities in the subsequent sections and partially supports the choice 
of the Sobolev space. 

First we assume that our function space Tt on £ is L 2 , i.e., 



H = U I / \(f)\ 2 d 3 x < oo}, 



and that the tentative domain T>(A) of the Laplacian A is C^°(E). Later we consider 
the case that 7i = H l instead of L 2 and see what is the difference. We do not claim 
that this analysis is new but we demonstrate this because we believe that this is the most 
illustrative explicit model in which the choice of the Hilbert space is highlighted. 

Separating the time variable t and the angular variables 6, ip we may write solutions 
in the form 

<f>im = e^f^YUO, <p) = e- ikt ^-Y lm (0, <p) (19) 

r 

with Yi m (9,(p) being the spherical harmonics. The reduced wave equation reads 

d 2 F t l(l + l)Fi 

-ri ~ 1 J + k 2 F l = 0, (20) 
ar z r z 

and the L 2 norm squared / dx 3 \(p\ 2 reduces to 

dr\F\ 2 (21) 



7 



up to an unimportant constant multiple. 

The behavior of the radial function near the origin is 



either T x ~ r l+1 (fr-r 1 ), (22) 
or ft~r-' (gt-r' 1 - 1 ). (23) 

All the Tis belong to the Hilbert space TC. The modes Q\ (I > 1) are not square integrable 
at r = and therefore are not normal modes. The mode Qi=q ~ const. (gi=o ~ cons t- ) is 
the only mode which requires further care. This mode is square integrable at r = 0. In the 



x 



case of £ ~ R 3 this mode does not belong to our Hilbert space because A(~) = — 47r5 3 ( 
is not in L 2 class. However, in the case of £ « R 3 — {0} this mode is allowed unless one 
further imposes a boundary condition at r = 0. However, the boundary condition to be 
imposed is not unique. Actually a boundary condition, 

aF' - F = 0, (24) 

is possible at the origin r = 0, where a is an arbitrary real parameter. In this sort of 
simple model one can immediately convince oneself that this is the most general boundary 
condition at the origin for the self-adjointness of the Laplacian operator but there is a 
systematic way to get the most general boundary condition, which is powerful for less 
simpler cases. We defer the demonstration of that method to the following subsections. 
Let us concentrate on the S-wave solutions (/ = 0). The most general S-wave solution 
which satisfies the above boundary condition is spanned by 

sin(fcr) ., . . , 

F k = — )—L + a ■ cos(Jfer) 25 
k 

with a being the constant in Eq. (|24]) . In this case we say that the self-adjoint extension 
is not unique so that the Laplacian with the initial domain T>(A) = C^°(R 3 — {0}) is not 
essentially self-adjoint and therefore the initial value problem is not well-posed. In the 
case of S ~ R 3 we have instead 

fi = (26) 

which contains no arbitrary parameter so that the Laplacian with the domain spanned 
by iVs of Eq. (|26|) is the only self-adjoint extension. | Therefore the Laplacian with the 
initial domain C^°(R 3 ) is essentially self-adjoint. 

This difference may be slightly more dramatic if we consider a spacetime M = S x R: 

-{0} fort>0, 
for t < 0. 




8 The extended domain spanned by F^s of Eq. ( pq ) is the if 1 closure ofC^°(R 3 ), and the corresponding 
self- adjoint extension is the Friedrichs extension S. 
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r=0 



t = 




■ o 



Figure 2: A conformal diagram of a spacetime M = £ x R with £t<o ~ R 3 and 
T, t>0 m R 3 — {0} so that a timelike singularity emerges at the center after t = 0. 

Namely, a timelike singularity emerges for t > as depicted in Fig. 0. In such a spacetime 
the normal modes do not match at t — unless a = so that the initial value problem is 
ill-posed. 

3.2 A systematic method of self-adjoint extension 

From the previous subsection we see that the problem of the function space L 2 (R 3 — {0}) 
for the field reduces to the problem of L 2 (0, oo) for the reduced radial wave function F. 
Let us study the solutions F± G L 2 (0, oo) of the equations: 

~^r = ±^±- (28) 

which are reduced from the equations (0) concentrated on the S-wave again. The solu- 



It is clear that the solutions Q± are not in L 2 (0, oo) class, while T± are. 

The prescription to find the most general boundary condition is to compose 



tions are 




(29) 



(30) 



F = F Q {r)+F + {r) + UF + {r) 



(31) 
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where F (r) G C^°(0, oo) satisfies the boundary condition F (0) = -Fq(O) = at the origin. 
U is the isometry of the space {J-"+} into the space {J 7 -} with respect to the L 2 (0, oo) 
norm, i.e., UJ 7 + (r) = e lOL T-{r). An elementary computation shows that 

F (°) = 6 V2 /om 

F(0) f + e ia { ' 

is a real number which we set equal to a -1 . This is what we alluded before. 

3.3 Sobolev space instead of L? 

Let us now change the Hilbert space from L 2 (R 3 — {0}) space to the Sobolev space 
if 1 (R 3 — {0}). We shall look for the solutions for which the integral 



drr 



o 



f±(r) 



2 r°° 



drr 



o 



'F±(r) 



(33) 



is convergent. However, we can see from Eq. (p5|) that the integral (|33| ) is divergent for 
T±. Therefore f± do not belong to the Sobolev space if 1 (R 3 — {0}) so that there remains 
no room to relax our boundary condition. That is, the Laplacian operator with the initial 
domain Cq°(R 3 — {0}) is essentially self-adjoint. 



Consider now the previous spacetime model: M=ExR (27) and the solutions of 
the field equation (]17D in the Sobolev space if x (E). It is now clear that the self-adjointly 
extended domain of the Laplacian agrees in both regions of the spacetime. Therefore 
the initial value problem is well-posed in the whole spacetime M. Actually the spherical 
wave propagates with no trace of the would-be singularity at the origin. Of course, this 
is because the spacetime is almost Minkowski. In a general wave-regular spacetime, the 
wave would be distorted and scattered by strong curvature there in a definite and unique 
way. 

It is physically assuring to see that the removed point is completely of no effect if the 
initial field configuration has a finite energy. This also supports that our choice of the 
Hilbert space is physically sensible. 



4 Wave probe in static spacetimes 
4.1 Spherically symmetric static spacetimes 

To illustrate the test of the essential self-adjointness of the operator A in Eq. (|T3|) in 
curved spacetimes, we first study the well-known spherically symmetric spacetimes. 
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In a general (n + 2)-dimensional spherically symmetric static spacetime, the metric is 
given by 

ds 2 = -V 2 dt 2 + V~ 2 dr 2 + R 2 dtt n . (34) 

Here we assume that V 2 is a positive function of r for < r < oo and is singular at 
r = so that the causal structure of the spacetime is as shown in Fig. [I]. Provided 
ip = f(r)Y(Cl), the equations (|l~6l) reduce to 

(V 2 R n Y c f 

where the prime denotes the derivative with respect to r and c is the angular momentum 
quantum number. The norm of / is given by 

,2 



Q_ 
2 



J dndrR n V~ 2 \f\ 2 + dfidrR n V 2 \f\ 2 , (36) 



where d\x is the volume element on the unit n-sphere (oE = dfidrV~ 1 R n ). 

For the essential self-adjointness of A, the norms of the solutions of the equations (|35|) 
should be divergent for each c and each sign of the imaginary term. We can easily verify 
that the norm 1 1/| | is divergent for c > if it is for c = 0, so we will examine the essential 
self-adjointness for the c = (S-wave) case. 

4.1.1 Negative mass Schwarzschild spacetime 

The 4-dimensional negative mass Schwarzschild metric is given by 

V 2 = 1 + , R = r, (M > , 37 

r 

and a timelike singularity is located at the center r = 0. Near the singularity, the equa- 



tions fl35D become 

/" + -/' = 0, (38) 
r 

since the other terms are less singular or even regular at r = 0. Then the two independent 
solutions / = const, and g = lnr are obtained. For the latter solution, the second term 
of the norm squared ([36D behaves as 

~ f drr 2 V 2 \g'\ 2 ~ lnr| -> oo. (39) 
Jo 

Thus the operator A on this spacetime is essentially self-adjoint, hence the spacetime is 
wave-regular. 

One might worry that our analysis does not involve the ±i part and only one of the two 
solutions for Eq. ([J5|) is verified not to be in the Hilbert space. One may be unhappy about 
the lack of intuition in the test being not completely convinced by the demonstration in 
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Sec. |3|. Here, we should remark that the solution / which well behaves near the singularity 
r = is divergent at infinity r = oo so that the Sobolev norm is divergent. The point 
is that if one of the independent solutions fails to well behave near the singularity there 
are no ways for the other solution to meet the condition at infinity because there is no 
available other independent solution to superpose. 

It is amusing to note that this also holds for the higher dimensional (n > 3) negative 
mass Schwarzschild spacetimes. Note also that if L 2 were chosen as the Hilbert space, the 
operator A would not be essentially self-adjoint in this case [(§]. 



4.1.2 Reissner-Nordstrom spacetime 

For the 4-dimensional over extreme Reissner-Nordstrom metric, 



2M ( Q 2 



V 2 = l + ^-, R = r,(Q 2 >M 2 ) } (40) 



where Q denotes the electric (magnetic) charge. Near the timelike singularity at r = 0, 
the equations (|35| ) for the S-wave (c = 0) become /" = (V 2 R 2 ~ Q 2 ) and / behaves as 
/ ~ r or a constant. Then, the norm squared (|36|) is finite. Thus, the classical singularity 
remains wave-singular. For c > modes, the norms are divergent. This implies that only 
the S-wave can fall into the singularity. 

Our analysis is basically local in time so that our wave approach to singularity probe 
can also work in the case that timelike singularities are hidden behind horizons by ex- 
tending the analysis in a straightforward way. 

Higher dimensional generalizations of the Reissner-Nordstrom solution are given by |], 



C D 2 

V 2 = l r + 7T , R = r, (41) 

with the parameters C and D being proportional to the mass and the charge, respectively. 
It is remarked that the singularities in the higher dimensional (n > 3, D 2 ^ 0) solutions 
are also wave-singular. 

To summarize the two examples above, the timelike singularity in the negative mass 
Schwarzschild spacetime is wave-regular while that of the over extreme Reissner-Nordstrom 
spacetime is not. The tendency that the over extreme Reissner-Nordstrom spacetime is 
more singular than the negative mass Schwarzschild spacetime sounds natural because 
the curvature is more divergent for the over extreme Reissner-Nordstrom spacetime. 

This reminds us of the well-known examples in quantum mechanics: the Coulomb 
potential and the r~ 2 potential problems in 3-dimensional space. The former is essentially 
self-adjoint and the latter is not if the r~ 2 potential is attractive and too strong ||. 
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One might guess from the two examples above that the quasi-local mass given (in 
4- dimensional case) by 

M local = -| {g^(d,R)(d u R) - 1} (42) 

would be finite as r — > in the wave-regular cases while it is infinite in wave-singular 
cases. However, this is not the case in other models as we shall see below. We shall 



propose an intuitive criterion for the wave- regularity in Subsec. [O. 
In general, if the metric functions of the metric ([34]) behave as 

R ~ r p , V 2 R n ~ r k , (43) 

near the singularity r = 0, the equations ( |3~5|) for S-wave (c = 0) become 

/" + = 0, (44) 
r 

under the condition np > k — 1, which holds for all our examples. Then, the solutions 
are / = const, and g = lnr (k = 1), or r 1 ^ k (k ^ 1). Since the norm squared fl3"6|) for the 
solution g = r (and similarly for g = lnr) is estimated as 

| \g\ | 2 ~ J drr~ k+2inp - k+1) + J drr- k ~ J drr~ k , (45) 

the singularity turns out to be wave- regular for the case k > 1. The locally flat example 
discussed in Sec. |3|is the case k = 2, the (n + 2)-dimensional negative mass Schwarzschild 
metric is the case k — 1 and the (n-|-2)-dimensional Reissner-Nordstrom metric = — ra+2. 

4.2 Other spacetimes 

We shall consider some less known but hopefully more physical solutions of the Einstein 
equations coupled to matter fields. The first two solutions below exhibit null naked 
singularities for some parameter regions. We shall remark on the null naked singularities 
from our point of view. 

4.2.1 The Wyman solution 

The Wyman solution is a static solution of the 4-dimensional Einstein equations coupled 
to a minimally coupled scalar field [O. The metric is given by 



m/n / \ m h 
r J \p + 2t]J 

= r 2(i_f!?) =p 1 - m ^(p + 2r ] ) l+m/ \ (46) 



r 



where n = \Jm 2 + a 2 with a scalar charge a, so m/n < 1 and p := r — 2n. A curvature 
singularity is located at p = 0. Since V 2 R 2 ~ p, this is the case of the metric functions 
with k — 1 and thus the spacetime is wave-regular. 
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4.2.2 Charged dilaton solution 

The 4-dimensional charged dilaton solution is given by |T0|j 



l-a z 



v = i-- r% («> 



r_ \ 1 + a 



2 a 



flf = r^l - - j , (48) 

with a 2 being a positive parameter in the model. Consider the extremal case r + = r_ . In 
the case a 2 > 1, the central singularity at p := r — r + = becomes timelike, while it is 
null for a 2 < 1. For any a 2 > 1, R 2 V 2 ~ p 2 this is the case of the metric functions fl43| ) 
with k = 2 and the singularity is wave-regular. The first term of the norm squared fl45|) 
diverges for a 2 < 3 so that we would reproduce the result in Ref. @ if we chose the L 2 
function space as our Hilbert space. 

4.2.3 String solution 

The 5-dimensional string solution given by 

ds 2 = V 2 (-dt 2 + dz 2 )+dr 2 + r 2 dn 2 2 , (49) 

has a curvature singularity at the center r = 0, which corresponds to a straight string. 
The operator A which appears in the wave equation □</> = — V~ 2 (d 2 +A)(j) = is expressed 
as 

a = _a; _ v* + + fj a r - ± } , (5i) 

where c is the angular momentum quantum number on the unit 2-sphere. By the separa- 
tion of the variables if) = f(r)e tkz Y(Q), the equations (A* =p = reduce to 

/" + _/' = 0, for S-wave, (52) 
r 

near r = 0. In this region, the solutions are / = const, and g = r~ 2 , and the norm 
squared for g 

\\g\\ 2 ~ j drr 2 r~ 4 + J drr 2 V 2 r- % (53) 
diverges. Hence the central singularity is wave-regular. 
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4.3 A simple criterion of wave- regularity for spherically sym- 
metric cases 



j dX 


df 




dX 



We may give an intuitive but not necessarily mathematically rigorous explanation of the 
wave-regularity and a simple criterion in what case the classical singularity becomes wave- 
regular. Take the example of static and spherically symmetric spacetimes. We can see 
that if we introduce a new radial coordinate X as 

the equations in the test of the essential self-adjointness look like 

_ c R^-2y2 f ± lR 2n f = 0> (55) 

and the essential part of the Sobolev norm is 

Af 2 

(56) 

Therefore the problem becomes similar to the quantum mechanics in a semi-infinite region 
(except the norm) if the variable X G (a, oo) for r G (0, oo) with a being a finite number. 
As is well known in that case the essential self-adjointness becomes non-trivial. On the 
other hand, if the variable X G (—00,00) for r G (0,oo), such a "half-space problem" 
would not appear. For the wave-regular case such as the negative mass Schwarzschild 
metric the range of X extends to —00 as the singularity r = is approached while it is 
finite for the over extreme Reissner-Nordstrom metric, which is the wave-singular case. 
Indeed this holds for all the cases given by the metric functions (f43|). This may suggest 
that the variable X in the wave mechanics plays a role similar to the affine parameter in 
the particle mechanics and that for a wave in a wave-regular spacetime the singularity is 
effectively infinitely far away. 

As a byproduct of the above observation we can see that if R 2 < 00 as r — ► and 
the singularity is null, i.e., / drV~ 2 —>■ 00 as r — > 0, then the singularity is wave-regular 
because / drV~ 2 R~ n — ► 00. This can be checked in the Wyman solution replacing a by 
io so that the parameter m/rj > 1 and therefore the singularity is null, though the scalar 
field becomes a ghost and the model becomes unphysical. In the charged dilaton model, 
the singularity becomes null for a 2 < 1 and is wave-regular. The wave-regularity of the 
null singularity was also asserted in Ref. Gj] but with different reasoning. 



5 Generalizations 
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5.1 Probes with other fields 



So far, the examples of timelike singularities have been probed with minimally coupled 
massless scalar fields. One can think of probing singularities with other fields such as 
spinor, vector, tensor fields, or metric perturbations. Our procedure of probing singularity 
is immediately generalized to each case by replacing the inner product (|T2|), hence the 
norm, in an appropriate way for the probing field. 

For example, when probing a timelike singularity with a massless scalar field coupled 
to the scalar curvature, we should adopt the inner product respecting the stress-tensor T^ v 
for the field (see Eq. (3.190) in Ref. lfL2|| ). The field equation in the case is (D + ^TZ)4> = 0, 
so the operator 

A = -V&VDi - iV 2 TZ (57) 

should be examined for the essential self-adjointness, where £ is a numerical factor and 
71 the scalar curvature. In the case of the Wyman solution, the equations (A* =p i)(f> = 
reduce near p = to 

p p l ZT) 2 

and the solutions are given by / = Ae l ^ ln p + Be~ l ^ ln p , (.4., B G C). After some calcula- 
tion, it is observed that the Sobolev norms for the solutions logarithmically diverge near 
p = 0. Thus, the singularity of the Wyman solution is also wave-regular when probed 
with the scalar field coupled to the scalar curvature. As is well known, in the confor- 
mally coupled scalar field case, that is £ = (d— 2)/4(d — 1) for any spacetime dimension 
d, the field equation is invariant under the conformal transformations of the metric and 
the field, g^ u (x) — ► g^ v (x) = C 2 (x)g^(x), <p — ¥ 4> — C^ 2 ~ d ^ 2 (j). Since T^„ transforms 
as T^ u = C 2 ~ d T^ u , the corresponding inner product is conformally invariant. Thus, the 
calculation is as simple as that in the static case when singularities in conformally static 
spacetimes are probed with conformally coupled scalar fields. This will be seen in the 
following subsections. 

Generalizations to the probes with spin 1/2 and 1 fields are similar. Explicit expres- 



f' + -J' + if = <>> (7=^), (5* 



sions of Tav for such fields are found, for example, in Ref. fl2fl 



5.2 Conformally static spacetimes 

In the previous section, we have examined the would-be naked singularities in static 
spacetimes, which can be regarded as regular if they are not detectable by waves, i.e., if 
the initial value problem for the field is well-posed. However, physically more interesting 
cases which can confront with the cosmic censorship are such dynamical spacetimes that 
naked singularities emerge after gravitational collapse. For example, to get insights into 
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the problem of the final fate of gravitational collapse, the Tolman-Bondi solution has been 
extensively studied by many people. It has been revealed that for some initial data shell- 
focusing naked singularities can be formed as a final product of spherical dust collapse JT^| . 
As the wave approach to fully dynamical cases is still far reaching for us at present, we 
will present a simple model to mimic such a dynamical problem on the basis of our study 
of the timelike singularity in conformally static spacetimes. 



5.2.1 Self-similar case 

Dynamical problems such as the formation of naked singularities can be made tractable by 
assuming self-similarity and the nature of naked singularities have been investigated |TJ] . 



Self-similar spacetimes also have attracted attention in connection with the critical be- 



havior in gravitational collapse and have been studied in detail by many authors [15]. 
Probing timelike singularities especially in self-similar spacetimes therefore is an interest- 
ing issue concerning the cosmic censorship. A technical advantage of self-similar metrics is 
that they can be written in the conformally static form so that we can straightforwardly 
apply our procedure developed for the static spacetime in the previous sections to the 
self-similar spacetimes when probing with conformally coupled scalar fields. 

For a massless scalar field, there is a non-static spherically symmetric solution discov- 
ered by Roberts |16 |. As one of the models of naked singularities in self-similar spacetimes, 



we will analyze the timelike singularity in the Roberts solution, whose metric can be writ- 
ten in the conformally static form, 

ds 2 = e 2ri ds 2 = e^{-d V 2 + dr 2 + R 2 (r)dn 2 }, 
R 2 {r) := l{l +p -(l- p )e-*}{<*-l), (59) 

where p is an integration constant and {d v )^ is the homothetic vector. For the value 
< p < 1, the curvature singularity located at r = becomes timelike and the global 
structure is identical to that of the negative mass Schwarzschild spacetime. By using this 
solution the problem of self-similar scalar field collapse has been discussed ]T6], |T7J. 

When probing the singularity with a conformally coupled scalar field, we can carry 
out the previous analysis with respect to the static metric ds 2 instead of ds 2 . Since near 
r = 0, R 2 ~ r, this is the case of the metric functions ( pE5| ) with k — 1. Thus, the 
norm squared (flop logarithmically diverges so that the singularity is wave-regular for the 
conformally coupled scalar field. 



5.2.2 Conformally flat spacetimes with emerging naked singularity 

More physical situations of gravitational collapse require that the spacetime contains 
a regular initial spacelike hypersurface on which the collapsing matter has a compact 
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support. 

We can construct a spacetime which has a regular initial hypersurface S io and a 
timelike singularity being formed in the future of S to . For example, let us consider a 
conformally flat metric 

ds 2 = C 2 {-dt 2 + dr 2 + r 2 d£l 2 } , (60) 
with a conformal factor which behaves near the center r = as 

r' (I > 0), 
r p + l (t<0), 

where p > —2 and C 2 (r = 0,t) — > sufficiently smoothly as t — > 0. Then, the spacetime 
has a timelike curvature singularity at the center r = for t > as depicted in Fig. |^. 
However, for sufficiently remote past, there are regular hypersurfaces S 4 . Therefore we 
can take an initial regular hypersurface E to at some to < and construct the Hilbert 
space on £ to . 

Since all the hypersurfaces S 4 are isomorphic to the initial hypersurface S <0 up to a 
conformal factor, the Hilbert spaces of a conformally coupled scalar field on E t are the 
same, even when H t intersect the central singularity. Therefore the spacetime is wave- 
regular for the conformally coupled scalar field. | 

Timelike singularities of the type examined above turn out to be also wave-regular 
when probed with the Maxwell field since it is conformally invariant and Maxwell's equa- 
tions are reduced to that of a massless Klein-Gordon field under a suitable gauge condition. 



5.3 Cylindrically symmetric case 

So far, we have studied spherically symmetric examples, whose central singularities are 
thus considered to be point-like. Here, as another example, we will probe a singularity of 
a cylindrically symmetric spacetime given by the metric, 

ds 2 = -p 2r71 dt 2 + p 2a2 dz 2 + p 2a3 dip 2 + dp 2 , (62) 

where the parameters cr, satisfy J2i a % — HiiPi) 2 — 1- This metric is known as the timelike 
Kasner solution and describes a cylindrical vacuum spacetime with a timelike curvature 
singularity along the line p = for the parameter values |crj| < 1. 

9 An everywhere smooth example can be constructed by choosing the conformal factor as 

C 2 = f U> x e-V*- 2 (fori>0), 
[U<x (e- 1 ^ 2 + e -!/t 2 ) ( for t < 0), 

where U> and U < , the analytic functions of r and t, should be chosen so that C 2 — > 1 as r — > oo. By 
taking U's appropriately, we have regular hypersurfaces at t < whose most of portions are flat. 
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On each hypersurface S t , the operator A in Eq. ( |T3| ) is written as 



(63) 



Since |<jj| < 1, the equations (p~6|) reduce to 




(64) 



near the singularity p = 0. In this region, the solution behaves like / = const, or g = In p 
and the norm squared becomes 



For the solution g = lnp, the second term logarithmically diverges, thus the singularity 
is wave-regular. 

For Oi = (0,0,1), the spacetime is flat and, if the angle coordinate <p has a deficit, 
the line p = becomes a cone singularity, which thus expresses a thin cosmic string in a 
locally flat spacetime. Also in this case the spacetime is wave-regular. 

6 Hair of wave- singular naked singularities 

In the previous sections, we have seen that most of the timelike "singularities" are wave- 
regular and that the singularity of the Reissner- Nordstrom spacetime is the only exception 
in our examples. In the wave-singular case, since probing waves feel the existence of the 
singularity in some sense, we naively expect that the waves scattered by the singularity 
will inform us of some feature of the singularity. Here, we shall discuss how to characterize 
wave-singular naked singularities. 

In the wave-singular case, the symmetric operator A in the wave equation has many 
different self-adjoint extensions. As discussed in Subsec. [27^, each self-adjoint extension 
corresponds to a different boundary condition at the singularity and accordingly describes 
a different time evolution of the wave for the same initial data. In other words, a wave- 
singular naked singularity has degrees of freedom for the possible choice of the time 
evolution operator. The degrees of freedom can be interpreted as the character or the 
hair of the wave-singular naked singularity. Since the self-adjoint extensions are in one- 
to-one correspondence with the set of partial isometries between the deficiency subspaces 
KL± of A (see Appendix for the definition), the set of isometries U describes the degrees 
of freedom. 

In general, when probing timelike singularities in static spacetimes with massless Klein- 
Gordon waves, we will find a positive real symmetric operator of the form A = —VD % VDi 




(65) 
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with domain C^°(E). Then the deficiency indices n+, n_ of A are equal and self-adjoint 
extensions can be made. If the naked singularity is wave-singular, then n + = n_ = N ^ 
and hence the partial isometry U is represented by an N x N unitary matrix U(N). Then, 
we say that the singularity has a U(N) hair. 

Let us consider a wave-singular spacetime which has an asymptotically flat region and 
a static region in a neighborhood of the central timelike singularity. More precisely, we 
will consider such a spacetime that the metric form is given by Eq. fl34|) and the metric 
functions behave as Eq. fl43|) with k < 1 near r = 0. Consider the wave probe with a 
massless Klein-Gordon wave, in particular the S-wave, in this spacetime. The solutions 
of the equations (A* =p i)ip± = near r = are given by 

f ± ~a ± r 1 - k + b ± , (66) 

where ij)± = f±(r)Y(Q) and a±, b± are constants. Both f± well behave near r = 0. On 
the other hand, in the asymptotic region, the equations (A* =f i)ip± =0 reduce to 

^ + ? /i + { ±i -^i)} /± = 0, (67) 

and the solutions for S-wave (I = 0) are exp {(1 =p i)r/\/2| /V and exp {(—1 =F *)^/v^2| 
Clearly the former solutions diverge in the asymptotic region and hence do not belong to 
the Hilbert space while the latter span the deficiency subspaces. Therefore the deficiency 
indices for I = mode are (1, 1) and the isometry is U(l). Thus, the wave-singular naked 
singularity of this spacetime has a U(l) hair. 

More detailed study will be given in the future work, in which we shall investigate 
what might occur in such a wave-singular spacetime in quantum field theory. 



7 Summary and Discussion 

We have studied the well-posedness of the initial value problem of the scalar wave equation 
in the case of static and conformally static spacetimes with timelike singularities choosing 
the Sobolev space as the natural Hilbert space. The physical idea behind the choice 
is that we can prepare an initial data only with finite energy. We have examined in 
detail the essential self-adjointness of the operator A in the wave equation defined in the 
Hilbert space in various models of spacetimes which contain timelike singularities in the 
conventional sense. In the spacetimes like the negative mass Schwarzschild spacetime the 
classical singularity becomes regular if probed with waves while more stronger classical 
singularities like the over extreme Reissner-Nordstrom spacetime remain singular. 
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We should comment that the wave-regularity of a spacetime does not guarantee that 
the spacetime is physically realizable. For example, a negative mass Schwarzschild space- 
time is wave-regular but allowance of negative mass solutions would make Minkowski 
spacetime unstable as pointed out by Horowitz and Myers (T^j. Probably there is a 
physics which rules out the negative mass Schwarzschild solution. In contrast, we may 
say that the over extreme Reissner- Nordstrom spacetime is unphysical on the ground that 
it is wave-singular. 

We have briefly touched upon the case that the timelike singularity emerges at some 
point in spacetime in rather artificial models which are the Minkowski spacetime with a 
single spatial point removed and the spacetime model which is conformal to Minkowski 
spacetime. However, timelike singularities of general spacetimes without any (conformal) 
timelike Killing vector will be more interesting from the view point of the cosmic cen- 
sorship especially in the case of gravitational collapse. In the case that spacetimes with 
"naked singularities" as illustrated by Fig. |3] are found to be wave-regular, such space- 
times can be said to be "globally hyperbolic" in the sense that the initial value problem 
is well-posed with the unique boundary value at the classical singularity. Clarke gave 
a sufficient condition for the well-posedness of the initial value problem for test fields. 
However, it turns out that the curve integrability condition for a naked singularity in 
Ref. is not satisfied for almost all the wave-regular cases discussed in the present work. 
We suspect that the condition is too demanding for the well-posedness of the initial value 
problem and his theorem can be further sharpened. 

The application of the present work to quantum field theory in curved spacetime is 
most interesting. The normal modes are solutions of the wave equation and an analogue 
of the Klein-Gordon inner product exists; i{(f, d t g) — (d t f, g)} with (/, g) being the inner 
product providing the Sobolev norm, which conserves if the spacetime has a (conformal) 
timelike Killing vector. 

We may carry out similar analyses to spinor, vector, and tensor fields. A natural 
question will be: Does the wave-regularity depend on with what fields we probe? At the 
moment what we can say is that in principle yes, i.e., it depends and it is nothing wrong 
from physical point of view. If initial value problems are well-posed for all fields, "naked 
singularities" are harmless and nothing to be afraid of. 
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Figure 3: A spacetime with emerging would-be naked singularities. 
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Appendix: Extensions of symmetric operators 

We briefly review some definitions of linear operators on a Hilbert space with an inner 
product (•,•). An operator on a Hilbert space H is a pair of a linear mapping A : H — > 7i 
and its domain of definition V(A). The pair (A,V(A)) is often abbreviated by A. If an 
operator A with V(A) densely defined in H satisfies 



then A is called symmetric. In the case, any vector v G Ti. can be approximated by vectors 
in T>(A) as close as possible. An operator A' is called an extension of A, if T>(A) C D(A') 
and A'^ = Aip, V0 G 22(A). Extensions of an operator A are obtained by the relaxation 
of the boundary condition on T>(A). Consider sequences {ip n } C T>(A) such that there 
exist limits lim n ^ oc -?/' n =: £ G 7i and lim^oo A-^ n =: ( <E H. If, for every such sequence, 
£ G £>(A) and A£ = (, then is said to be closed. If a non-closed operator A 

has a closed extension it is called closable. Every closable operator has a smallest closed 
extension, which is called its closure. Consider a symmetric operator (A,V(A)). Define 
V(A*) to be the set of all G H for which there exists x £ W such that 



(0,^) = (A^), V0,^GP(A), 



(68) 



(0,^) = (x,^), V^6D(4 



(69) 
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Then, since T>(A) is dense, x is uniquely determined by G V(A*) and Eq. fl69|). An 
operator (A* ,V(A*)) defined by A*cf) = x for every G T>(A*) is called the adjoint of 
(A, T>(v4*) may be larger than T>(A), in which case A* is a proper extension of A. 

If (A*,V(A*)) = (A,V(A)), an operator (A,V(A)) is said to be self-adjoint. 

Now let us see an example of extensions of symmetric operators to self-adjoint ones. 
Take Ti = L 2 (0, 1) and consider an operator T: 

Ti, = -»*^>, (70) 
V{T) = {i-\i,(0) = 4,(l)=0,i-€ACl0,l]}, (71) 

where AC[0, 1] expresses the set of absolutely continuous functions on [0, 1] whose deriva- 
tives are in L 2 (0, 1). It can be verified that T is symmetric. For 0(x) = exp(ikx) G T>(T*), 
T*4> — X — kexp(ikx) G Ti is not an element of T>(T), thus V(T*) D V(T); T is not 
self- adjoint. 

Next, consider an operator T a with the same action as T in T>{T) and with the domain 
V(T a ) := {> | V(0) = e fa V(l), V> e AC[0, 1]}, (72) 

where a is a real number. Clearly, this is an extension of T. For <p G T>(T a *), there exists 
X G 7i such that 

WGP(T a ), (<P,T a iJ>) = (x,il>), T a * ( f ) = x. (73) 

Namely, 

cix0* [ -i^ ) = ( l dx (T Q *0)* ^. (74) 
\ dx ) Jo 

On the other hand, by partial integration it is obtained 



<£z (T a *0)* V- (75) 



Thus, <p*(l)ij}(l) — 0*(O) , 0(O) = 0. From the boundary condition for ip, it follows that 
0*(1) = e ia (j)*(0) = (e- ia (f)(0))% hence 0(0) = e iQ 0(l). Therefore V(T a *) = V(T a ); T a 
is self-adjoint. Since a is arbitrary, it turns out that T has infinitely many different 
self-adjoint extensions. 

In general, for a closed symmetric operator A, the closed symmetric extensions can 
be carried out in more systematic way. Consider solutions of A*<f>± = ±«0± and the sets 
K,± of the solutions 0±, respectively, which are called the deficiency subspaces of A. The 
pair of numbers (n + ,n_) := (dim/C + , dim/CL) is called the deficiency indices of A. If 
A is a closed symmetric operator with the deficiency indices n + = n_, then A has self- 
adjoint extensions, and if n + = n_ = 0, A is self-adjoint. Let U be the partial isometries 
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/C+ — > /C_. Then, the self-adjoint extensions can be obtained by taking the domains 



as 



V(A E ) := {0o + 0+ + | 0o G V(A), 0+ G £+}. (76) 
In the example above, the normalized solutions 0± of the equations T*<f>± = ±«0± are 

4> ± (x) = ] Ze (77) 



The deficiency subspaces of T are /C± = {P4>±\ P G C} and thus the deficiency indices 
are (1,1). Then, the partial isometries are taken as U : /C+ — > /C_ : 0+ i— ► 70- where 
|7| = 1. The symmetric extension with respect to 7 G £7 is given by T 7 = —id/dx with 
the domain 

P(T T ) := {0 O + /5(0 + + 7 0_) I 0o G P(T), /5 G C}. (78) 
It turns out that the phase factor e ta is given by 

&ia = m = l±7e_ (7Q) 
0(0) e + 7 V ; 

If the closure of a closable symmetric operator A is self-adjoint, A is called essentially 
self-adjoint. In this case, A has a unique self-adjoint extension. The basic criterion for 
essential self-adjointness is to verify that its deficiency indices are both zero, namely, the 
solutions 0-t of the equations (A* z fi)4>± = are not in the considering Hilbert space. An 
example of essentially self-adjoint operators is the Laplacian operator on L 2 (R n ) with the 
domain C JO (R ri ): a set of smooth functions with compact support. 

Detailed studies of extensions of symmetric operators and further examples can be 
found in the text book HI. 
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